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Abstract 

The pion mass difference generates a pronounced cusp m K ^ 3-k decays. As has 
recently been pointed out by Cabibbo and Isidori, an accurate measurement of the 

QQ . cusp may allow one to pin down the S-wave vrvr scattering lengths to high precision. 

^^ I Here, we present and illustrate an effective field theory framework that allows one 

to determine the structure of this cusp in a straightforward manner. The strictures 

\^ . imposed by analyticity and unitarity are respected automatically. 

o 
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1. The S-wave vrvr scattering lengths 00,^2 have been predicted with percent 
5^ \ level accuracy some time ago [1,2], and first steps for an experimental veri- 

fication of this prediction have been performed in Ref. [3,4]. Recently, it has 
been pointed out by Cabibbo and Isidori [5,6] that isospin violating effects 
generate a pronounced cusp in if — >■ Svr decays whose experimental investiga- 
tion may allow one to determine the combination oq — 02 with high precision. 
A first analysis of data based on this proposal has appeared [7]. (The strong 
impact of the unitarity cusp on vr'^vr'^ scattering close to threshold was already 
mentioned in [8].) In order for this program to be carried through successfully, 
one needs to determine the structure of the cusp with a precision that matches 
the experimental accuracy. In view of the large amount of data available [7], 
this is a considerable task. In the present letter, we present a method which - 
we believe - has the potential to achieve this goal. 
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In [5,6], the structure of the singularity at the cusp is investigated using 
unitarity, analyticity and cluster decomposition properties of the S'-matrix. In 
addition, an approximation scheme is used, which consists in expanding the 
decay amplitude in powers of tttt scattering lengths. The latest work [6] retains 
effects up to order (scattering lengths)^ and omits explicit electromagnetic 
effects. Here, we present a Lagrangian framework, which automatically satisfies 
unitarity and analyticity constraints and, in addition, allows one to include 
electromagnetic contributions in a standard manner. Specifically, we use a non- 
relativistic framework that has already proven to be useful in the description 
of bound states [9-25]. In this framework - in contrast to relativistic field 
theory - an expansion in powers of scattering lengths emerges automatically 
from the loop expansion. Moreover, it is a scheme that provides a proper power 
counting. [The low-energy expansion proposed here is closely related to early 
work performed in the sixties by many authors (for a review see [26]), who used 
S'-matrix methods to investigate the production of particles - in particular also 
in i^ ^ Stt decays - in the threshold region. The method presented here may 
be considered an effective field theory realisation of these approaches.] 

The strategy that we follow in this letter is the following. First we write down 
the most general non-relativistic Lagrangian relevant for this decay and deter- 
mine all four-pion couplings therein through a matching procedure in terms of 
the threshold parameters of tttt scattering. In the next step, we evaluate the 
if — > Stt decay amplitudes to two loops. This results in an explicit represen- 
tation of the S'-matrix elements which is valid in the whole decay region and 
slightly beyond. We propose to analyse the experimental data with the use 
of this representation. We plan to include real and virtual photon corrections 
at a later stage. On the other hand, we do keep the pion and kaon masses at 
their physical values, which is a fully consistent procedure in this framework. 
This guarantees that the various branch points and cusps occur at the proper 
place in the Mandelstam plane. 

We display the results without a detailed derivation, which will be provided 
in a forthcoming publication [27]. 

2. We consider the neutral and charged decay modes K'^{px) -^ 
7r°(pi)7r°(p2)'?'"^(P3) and K'^{pk) -^ 'n'^ {Pi)t^~^ {P2)'^~ (Ps) ■ The kinematical 
variables are defined as usual: Si = {Pk — PiY with pf = Mf , i = 1,2,3, 
where M^+ = M^^ and M^^o denote the masses of the charged and neutral pi- 
ous, respectively, and A^r = M^ — M^o 7^ 0. In the centre-of-mass (CM) frame 
Pk = {Mk, 0), with Mx the charged kaon mass. 
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where \{x,y,z) = x + y + z — 2xy — 2xz — 2yz is the triangle function. 



Below we also use the velocities Vjk and kinetic energies Tj, 



A(g.,M|,M|) , 

^ife(si) = i , Ti=p1-Mi . (2) 



3. A non-relativistic approach to describe decays i^ — i> Stt can be justified, 
if the typical kinetic energies of the decay products are much smaller than 
the masses. This can be achieved by considering a world where the strange 
quark mass is taken to be smaller than its actual value. Then, a consistent 
counting scheme arises, if one introduces a formal parameter e and counts Tj 
as a term of order e^, the pion momenta as order e, whereas the pion and kaon 
masses are counted as 0(1). From X]j Ti = Mk — J2i ^ii oiie concludes that the 
difference Mk — I]i ^i is then a quantity of order e^ as well. The pion mass 
difference A^ is also counted as O(e^). The effective field theory framework, 
which we construct below, enables one to obtain a systematic expansion of 
the amplitudes in e. For sufficiently small rris, the expansion in e is expected 
to work very well. 

Together with e, our theory has another expansion parameter, namely a char- 
acteristic size of the tttt threshold parameters, which we denote generically as 
a. In particular, the amplitudes in the non-relativistic framework are given 
in form of an expansion in several low-energy couplings Ci,Di, which can be 
expressed in terms of the threshold parameters of the relativistic tttt scattering 
amplitude. We expect the expansion in a to converge rapidly because of the 
smallness of the scattering lengths. These two expansions are correlated [27]: 
because one-loop integrals are of order e, adding a pion loop generated by a 
four-pion vertex increases both the order in a and in e by one. A consistent 
power counting is achieved: to a given order in a and in e, a well-defined finite 
number of diagrams contribute. 

Increasing now rris to its physical value again, convergence in the e-expansion 
is not a priori evident, because Ti/Mi can become as large as 0.4, and the 
corresponding maximal momentum |p| is then not much smaller than the 
pion mass. However, let us note the non-relativistic framework is only used 
to correctly reproduce the non-analytic behaviour of the decay amplitudes 
in the kinematical variables si, 82,83, and to thus provide a parametrisation 
consistent with unitarity and analyticity - a trivial polynomial part in the am- 
plitudes can be removed by a redefinition of the couplings in the Lagrangian. 
In addition, from the analysis of the experimental data one knows [6] that in 
the whole physical region the real part of the decay amplitude can be well 
approximated by a polynomial in 81,82, 83 with a maximum degree 2. We in- 
terpret this fact as an experimental indication for a good convergence of the 
e-expansion for the quantities one is interested in. 

We now proceed with the construction of the non-relativistic Lagrangian 
framework. In the decay amphtudes, we shall restrict ourselves to terms up to 



and including 0(e^, ae^, a^e^). 

4. It is convenient to formulate the non-relativistic approach in a manner 
that describes the two-particle subsystems in a manifestly covariant way [27]. 
We start with the tttt interaction and consider the following five channels in 
TT^TT^ ^ n^TT^: {ah- cd) = (1) (00; 00), (2) (+0; +0), (3) (+-; 00), (4) (+-; +-), 
(5) (++; ++)• The Lagrangian takes the form 

C^^ = 2Y,c!>lw±{tdt - W±)<^± + 2<^lWo{tdt - Wo)% + E A, (3) 
± 1=1 



where <l>j is the non-relativistic pion field operator, W± = JM^ — A, Wo 
M^o — A, with A the Laplacian, and 



/:,=X,a(<l>l<l>l<l>a<l>b + h.c) + X,Di\^{W,^lWd^\^a^b + ^l^\Wa^aWb^b 
+ V^lV^\^a^b + ^l^\V^a^<^b - h<^l<^\^a^,) + h.c] + . . . , (4) 



with hi = 2M^o, h2 = 2M,M,o, /i3 = 3M2 - M^o, h^ = h^ = 2Ml The 
ellipsis stands for terms of order e^ as well as for P-wave contributions, which 
occur at order e^ in the tttt amplitude. They do not enter the i^ — > Svr matrix 
elements at the order of accuracy considered, so we omit them here. The 
low-energy constants Ci, Di are matched to the physical scattering lengths 
below. To simplify the resulting expressions, we have furthermore introduced 
the scaling xi = x^ = 1/4, X2 = x^ = X4^ = 1. Finally, note that we omit 
local 6-pion couplings as well. Their contribution to the K ^ Sn amplitude is 
purely imaginary in the non-relativistic framework, and of order e^, see also 
[6,27]. 
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The pion propagator is given by 

(27r)4 2w,(p)(u;,(p)-po-?0)' 

(5) 
and u'-i-(p) = J Ml + p^, wo(p) = \l M^o + p^. The loops are evaluated by 
using the following prescription: one expands the square root in a series, cal- 
culates the emerging integrals using dimensional regularisation, and sums up 
the series. In this manner one ensures that the power counting in the non- 
relativistic theory is not destroyed by the loop corrections. 

All loops in TTTT scattering can be expressed through the basic integral 

^ _ f d^l 1 1 

■^'"'^ > ~ J Ji2^ 2wa{l)2w,{P - I) K(l)-/o)MP-l)-Po + /o) ' ^^' 
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which is a quantity of order e. 





(7) 



5. The couphngs Ci,Di can be expressed in terms of the threshold parame- 
ters of the underlying relativistic theory. In the isospin symmetry limit, the 
expansion of the relativistic amplitude reads 



Ref,(5,t) = A,{l + ^(s-4M^)} 



The ellipsis stands for higher orders in e, as well as for P-wave contributions. 
The bar indicates the isospin symmetric limit, at M^ = 139.57 MeV. In terms 
of the standard dimensionless scattering lengths Oq and a2, one has 

3Ai = N{ao + 2a2) , 2^ = Na2 , SA^ = N{a2 - ao) , 

6^4 = N{2ao + as) , A = Na2 ; iV = 327r , (9) 

with oq — ^2 = 0.265 ± 0.004 [1, 2]. The products Ajfj denote effective ranges. 
Still in the isospin symmetry limit, the couplings Cj, Di are related to these 
threshold parameters in the following manner, 

2a = Ai , SM^Di = A,n ■ (10) 

Taking isospin breaking into account, one finds at leading order in chiral per- 
turbation theory [28] 

2Ci,2,5 = A,2,5(l - V), 2C3 = ^3(1 + V3), 2C4 = ^4(1 + r^), (11) 

where 77 = A^/M^ = 6.5 x 10^^. Isospin breaking corrections in D^ do not 
contribute at the order considered here. 

6. It remains to display the i^ ^ 37r Lagrangian, 



Ck = 2K^Wk (idt - Wk]K 



+ i Go [k^<^+<^1 + h.c) + ^ Gi [k\W+ - M^)^+<^1 + h.c) (12) 

+ ^ i/o (i^^$_$+ + h.c) + ^ i/i (k\W^ - M^)<^_<^1 + h.c) + . . . , 
where K denotes the non-relativistic field for the K^ meson, Wk = 



M]^ — A, and the ellipsis stands for the higher-order terms in e. Note that 
all couplings Gj, H^ are assumed to be real. Their contribution to the decay 
matrix elements at tree level is provided below, in the amplitudes Ajsi c- 



The complete Lagrangian of the theory is Ck + Ct,.„. The tree-level expressions 
for the amplitudes, generated by Ck-, are modified by final state interactions 
of the pions, generated by loops evaluated with C,^t^. We use the notation 



-M00+ = M'^^^ + A^^-^°°P + A^ 2-loops ^ j^+ ^ ^0^0^+] ^ 

for the decay amplitudes and the Condon-Shortley phase convention for the 
pions. Our amplitudes are normalised such that the decay rates are given by 

''^-2]^(^'^''*'"(^'-^''i-^i\n^. (14) 

7. We now display the tree and one-loop results and modify the notation for 
the couplings Cj, Di in order to make the formulae more transparent: 

(Ci, C2, C3, C4, C5) = (Coo, C+o, Cx, C^-, C'+_|_) , (15) 

and analogously for the Di. We find 

Mn °°^ = Bni (S3) ^+- (S3) + ^ATa (S3) ^00 (S3) 

+ {5Af3(si)J+o(si) + (Si ^ S2)} , 

-M?^" = Ac;(s3), 

M^-^°°P=fici(s3)J++(s3) 

+ {5c2(si) J+-(si) + 5c;3(si) Joo(si) + (si ^ S2)} , (16) 

where ^ 



P?+P2 



7r J f ) 



^^(s3) = Go + Gi(p°-M.; 

5^i(s3) =2(a + /^.(S3 - sMh^ + i/i(^^^^ - M. 
Bn2{s^) = (Coo + I^oo(s3 - soo)) |Go + Gi(p° - M^)| , (17) 

E^3(si) = 2(C7+o + DUsi - s+o)) |go + Gi (^^^ (l + ^) -M.)| , 
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To render the formulae more compact, we keep some terms in Eqs. (17), (18) that 



contribute at order ae' . 




(A) 



(B) 



Q -1 , M , Q -k , M 



Fig. 1. Two topologically distinct non-relativistic two-loop graphs describing the 
final-state tttt rescattering in the decay K ^ 3tt, with Q^ = (^i + 92)^- 

and 



Bc2iSi) 

Bcsisi) 



C++ + D++(s3 - s++))\Ho + hM - M, 



2(C+_ + D+_{s, - s+_)){Ho + iJif^i±^ - M. 



a + D,{s, - s.))|g'o + G'i(p° - M^ 



In the above expressions, Sj denotes the physical threshold in the i channel 
and p° are given by Eq. (1). Note that, according to this equation, the masses 
in the relation of p° to Si differ in the neutral and charged channels. 

8. There are two topologically distinct two-loop graphs that describe pion- 
pion rescattering in the final state, see Fig. 1. At the order of accuracy we 
are working, it is sufficient to consider the case of non-derivative couplings. In 
this case, the contributions of both diagrams depend only on the variable s, 
where 

Q^ = (qi + q2r , Q' = s. (19) 

The diagram in Fig. IB, apart from a factor containing coupling constants, is 
given by a product of two one-loop diagrams which were already calculated in 
Eq. (6). The non-trivial contribution from Fig. lA is proportional to 



Mis) 



d^^l d'^k 



1 
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2wa{l + k) Wa{\ + k) - Mx + /o + A;0 2wi>{l) Wb{\) - P 
111 1 



X 



2wc(k) ti;,(k) - fco 2wrf(Q - k) Wd(Q - k) - QO + /^o ' 



(20) 



A detailed discussion of this integral will be provided in [27] . Here we simply 
note that the most general representation for diagram Fig. lA can be written 
in the form 

Mis) = f(m,, M,, M„ M,; s) + P(s) Jed(s) + P'{s) , (21) 



where F is ultraviolet-finite and contains the full non-analytic behaviour of 
the two-loop diagram in the low-energy domain. Further, J^d is the one-loop 
function Eq. (7), and P, P' are real polynomials. We have suppressed the 
dependence of F on the mass Mj generated by Q^, see below. 

In the following, we use a simplified form of F, where part of its imaginary 
part is dropped - this omission affects the decay width at order a^ only and 
is therefore of no relevance here. We use the integral representation [27] 



F[Ma,Mh,M,,Md;s) = 
where 



y^- ^K 1^^ {Ma-M,f\^'^ 



2v^V So ; (2(Mi + M2) - (M, + M,)2 - So 

2/1 , X\2\l/2 



1/2 



So = M'+ Mi - 2Mk{M' + 



g{y, s) = -(1 - y)ql - y^' + ''\,"'^X2,\^' " ^0 



2 , y{ l - y)Cl\l + 5f 
4(1 + |/Q7s) 



Qo = ^ , s = {Mc + Md) 



^,_\{Ml, Ml{M, + M,f) Ml - Ml 

AMI ' " s 



^. ^ ^^^,^^,^^_^_^^ ^ ^^ c d _ (23) 



Approaching threshold from above, we find 

go 

1287r2(M, + Mrf) So 



F{M.. M. M.. M. .) ^ - ,,,,.,:^,^/'''^°-5''^' + Oi.i) . (24) 



where So denotes So ^-t ^q = 0. This last relation shows that F is of order e^. 

9. Our prescription for the representation of the decay amplitudes at 0{a^) 
is as follows: we evaluate the contributions from all the graphs displayed in 
Fig. 2 and Fig. 3. Further, in the graphs of the type Fig. lA, we retain only 
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Fig. 2. Two-loop graphs contributing to the decay K^ -^ 7r'^7r'^7r+ in the 
non-relativistic effective theory. The graphs obtained by a permutation of identical 
particles in the final state are not shown. 

the non-analytic piece F, whereas the polynomials P, P' are included in the 
tree- level couplings Gi^Hi. This choice of a particular representation of F is 
equivalent to a renormalisation prescription. 

With this convention, we find for the amplitudes at order a?e^ 

_^ 2-loops ^ j^A^^^^^ ^^^ ^^^ ^ j^B^^^^ ^^^ s,); I = N,C , (25) 



where 



M^ = 4HoC+^C,F+iM^, M^, M^, M^; S3 



+ 2HoC++C^F+[m^, M^, M^, M^; sgj 
+ AG0C00C+0F+ (m^, M^o , M^o , M^o ; S3 ] 

4HoC^C+oFo{m^, M,, M,, M^o; Si 

+ AGoC%Fo[m^o, M^, M^o, M^; si] 

+ 2GoCooC+oFo(m^o,M^o,M^, M^o; sA + (si ^ S2) L 



A<« =4ffoaC+-4_(s3) + GoC^,Ji,{ss) 

+ 2 GqC^, + HqCxCoo J+_(s3)Joo(s3) 



(26) 



+ 
+ 





Fig. 3. Two-loop graphs contributing to the decay K~^ -^ Tr^ir^ir^ in the 
non-relativistic effective theory. The graphs obtained by a permutation of identical 
particles in the final state are not shown. 



M^ = 2GoaC++F_ M^o, M^o, M,, M^; S3 



+ i^4HoCl_F+ (m., M., M^, M.; s,] 

+ 2G'oaC+-F+(M,o, M,o, M,, M^; Si^ 
+ 2/7oC++C+_F+ (m^, M^, M., M^; s^ 

+ AGoC+oC^F+ (m,, M,o, M^o, M^o; si) + (si ^ S2; 



'■c 



iO'^+H 



IS3J 



J+-{si)Jooisi) 



+ 2 



HqC^ + G'oC2;C_|__ 



[Sl ^ S2, 



(27) 



Here, Fj(. ..;s) stands for the integral F(. ..;s), evaluated at Q^ = 
A(Mi, M2„ s)/4Mi, with ^ = ±, 0. 

10. The decay amplitudes depend on the four real i^ — > Stt coupling con- 
stants Hi, Gi and on the threshold parameters for tttt scattering. Combining 
the tree and one-loop result Eqs. (16)-(18) with the two-loop contributions 
Eqs. (25)-(27), we obtain the neutral and charged decay amplitudes up to and 
including terms of order e^,ae^ and a^e^, expressed in terms of the one- and 
two- loop integrals J and F displayed in Eqs. (7) and (22), respectively. This 
representation is valid in the whole decay region, and is the main result of this 
article. 
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The decay amplitude K^ -^ -k^-k^-k^ obeys what we refer to as the threshold 
theorem: the coefficient of the leading non-analytic piece, which is propor- 
tional to f_|__(s3), is given by a product of two factors, the decay amplitude 
K^ -^ TT+TT+TT^ and the scattering amplitude tt+tt^ -^ vr^vr^, both evaluated 
at threshold. Thus, the heuristic argumentation of Ref. [5], which serves as 
a cornerstone of the whole method, is conffimed in the effective field the- 
ory framework. This threshold theorem has its analogue in hadronic atoms, 
viz., in the modification of energy levels and decay widths through hadronic 
interactions [15-25]. Of course, aside from the determination of the leading 

term in v^ , our approach also allows a systematic evaluation of higher-order 

contributions v\_,v\_ .... 

11. We now compare the content of this letter with the recent work of Cabibbo 
and Isidori [6] (CI), who have proposed an alternative representation for the 
K ^ ^i-K decay amplitudes. Assuming certain analytic properties of the decay 
amplitudes, CI derive a representation of the amplitudes up to and includ- 
ing terms of order a^, using analyticity, unitarity, and cluster decomposition 
properties of the S matrix. 

Conceptual aspects of our methods have already been compared in the intro- 
duction. Here we add that, ffistly, we do provide an explicit representation 
of the decay amplitudes that is valid in the whole decay region, including all 
powers of velocities generated by the graphs considered. Secondly, we note 
that our explicit two-loop calculation conffims the expected [6] analytic prop- 
erties of the amplitudes in the vicinity of the tt+tt^ threshold. On the other 
hand, away from this threshold, but still in the physical decay region, branch 
points develop [27], contrary to the expectations spelled out in [6]. 

After these general remarks we compare the amplitudes in more detail at one 
and two loops. In the actual calculations in Ref. [6], an approximation has 
been used: the angular integrals have been replaced by averages, where the 
integrand is evaluated at a certain value of cos^. As CI note, this approxima- 
tion is exact, if the integrand is at most a linear function in cos^. Since this is 
true at 0{ae^), one expects that our results at this order algebraically agree 
with CI. We have checked that this is indeed the case (up to a few typos [29]). 

At two loops, our results are algebraically different from those of Ref. [6] . The 
reason can be traced back to the following. At two loops, the discontinuity 
cannot, in general, be obtained from an integration over cos^ without further 
ado: the path of integration has to be deformed into the complex plane - it 
does not simply run from cos^ = —1 to cos^ = 1 along the real axis [26]. In 
fact, near the pseudothreshold S3 = {Mk — Mj,)^, the deformed path runs to 
infinity, thus generating an infinity in the discontinuity there. In the case where 
all internal masses are equal to the charged pion mass, integrating cos 6 along 
the real axis does generate the correct discontinuity up to S3 = {Mj^ — M^)/2. 
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In the case where two neutral pions are running in the inner loop, the situation 
is more complicated, because an anomalous threshold develops in the lower 
half plane. Still, integrating along the path mentioned generates the correct 
discontinuity near threshold. [These difficulties do not arise in our approach 
- the function F in Eq. (22) is smooth on the upper rim of the cut, while its 
discontinuity develops the singularities mentioned, at the positions predicted 
by the Landau equations [27,30,31].] 

A replacement of the angular integrals by an average [6] , where the integrand is 
evaluated at a certain value of cos^, can therefore be a reasonable approxima- 
tion only in the vicinity of the threshold S3 = 4M^. Indeed, our result agrees 
(up to a typo [29]) with Ref. [6] to the order considered here, at threshold. 
Away from threshold, the expressions differ. 

Finally, we shortly comment on the (revised) article by Gamiz et al. [32], which 
appeared only very recently. It is the aim of that article to provide an error 
analysis of the procedure proposed in [5, 6]. The authors investigate the process 
ii' ^ 37r in the framework of chiral perturbation theory, and approximate 
two-loop graphs by retaining their discontinuity only - an approach which 
generates fake singularities in the transition amplitude, as just mentioned. 
An ad hoc prescription is invoked in [32] to avoid these singularities when 
investigating the cusp, see e.g. their comment after Eq. (4.20). 

12. In summary, we have investigated K ^ "in decays within a non-relativistic 
effective Lagrangian framework. The amplitudes are calculated in a systematic 
double expansion in the kinetic energies of the decay products (which we count 
as terms of order e^), and in the threshold parameters of elastic tttt scattering 
(which are generically denoted by a). We provide an explicit representation 
of the amplitudes at order e^, ae^, a^e^ - valid in the whole decay region - in 
terms of the (real) K -^ ^iti coupling constants Gj, Hi and of the threshold 
parameters a. 

Our amplitudes agree with the ones of Cabibbo and Isidori [6], up to and 
including terms of order at^ . On the other hand, at order a^, they differ away 
from threshold. We propose to repeat the analysis of the experimental data 
of the NA48/2 collaboration [7] with our representation of the amplitudes, for 
several reasons: i) In view of the aimed precision, one ought to examine the 
importance of the mentioned differences in the determination of vrvr scattering 
lengths, ii) It would be useful to extend the ffi to the full decay region, and 
to the charged decay modes K^ -^ tt+tt+tt^ as well, in order to determine a 
maximal set of tttt threshold parameters, iii) As was already pointed out in 
[6], cusps also occur at the border of the Dalitz plot. Investigating data in 
those regions may allow one to determine different combinations of scattering 
lengths. 
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It remains to investigate the importance of higher orders in the low-energy 
expansion, and to apply radiative corrections, which can be evaluated in the 
field-theoretical framework used here in a standard manner. The effects gener- 
ated by the tt+tt^ bound state at the tt+tt" threshold can also be investigated 
within the same approach (see, e.g., [15-23]). We plan to include these effects 
in forthcoming publications. 
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